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^v^j 1 Abstract. This is a study of the Euler equations for free surface water waves in the case of varying 

bathymetry, considering the problem in the shallow water scaling regime. In the case of rapidly varying 
'^j i periodic bottom boundaries this is a problem of homogenization theory. In this setting we derive a new 

model system of equations, consisting of the classical shallow water equations coupled with nonlocal evolution 
equations for a periodic corrector term. We also exhibit a new resonance phenomenon between surface waves 
' and a periodic bottom. This resonance, which gives rise to secular growth of surface wave patterns, can be 

i viewed as a nonlinear generalization of the classical Bragg resonance. We justify the derivation of our model 

with a rigorous mathematical analysis of the scaling limit and the resulting error terms. The principal 
issue is that the shallow water limit and the homogenization process must be performed simultaneously. 
0-H ' Our model equations and the error analysis are valid for both the two- and the three-dimensional physical 

problems. 

1. Introduction 

Studies of the Euler equations for free surface water waves are important to understanding the dynamics 
of ocean waves. The case of an idealized flat bottom and the resulting model equations has been widely 
studied for many years. The more realistic situation of varying bathymetry is less well known, despite its 
iy) | fundamental importance to studies of ocean wave dynamics in coastal regions, and there is not a complete 

consensus as to the appropriate model equations. In the case of topography there are many asymptotic 
| scaling regimes of interest, including long-wave of modulational hypotheses for the evolution of the free 

■ surface, and short scale and/or long scale variations in the variable bottom fluid boundary. 

In this paper we address the evolution of waves in the shallow water regime, for which we investigate 
the effect of the roughness of the bottom topography. The simplest situation is where the bottom varies 
periodically and rapidly with respect to the typical surface wavelength, a regime which can be described in 
the context of homogenization theory. Ideally, wave motion in this regime of rapid periodic bottom variations 
is described in terms of a long wave effective component, which is then adjusted by a smaller multi-scale 
corrector at the next order of approximation. In terms of the initial value problem for this regime, initial 
configurations consisting of large scale data with a multi-scale corrector term are expected to give rise to 
solutions with the same character, up to a smaller error term. In this paper we derive a system of model 
equations for such multi-scale approximate solutions. While other authors have looked at similar situations, 
as far as we know this system is new, consisting of a version of the shallow water equations for a mean field or 
effective components of the surface elevation and the fluid velocity, which then drive a nonlocal system of two 
additional equations for the evolution of a more rapidly oscillating corrector term. Because of the number 
of other models that have been proposed to describe this setting, we justify the derivation of our system 
with a rigorous analysis, giving error estimates for our approximate solutions. In cases in which there is a 
resonance between the effective velocity and the periodic bottom, the solution of the corrector equation can 
exhibit secular growth at a linear rate. This phenomenon can be viewed as a nonlinear generalization of the 
classical Bragg resonance between the bottom topography and the free surface. This is a local phenomenon, 
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which may occur when the local Froude number is subcritical. In the absence of resonances, our analysis is 
valid over time intervals of existence of the effective component. 

The literature on models of free surface water waves over a variable depth is extensive, including the 
paper of Miles [21] on its Hamiltonian formulation, and that of Wu [29] on models which are valid in long 
wave scaling regimes. The paper of Rosales & Papanicolaou [27] studies the long wave regime in which the 
bottom is rapidly varying, in the sense that the typical wavelength of surface waves is taken to be much 
longer than typical lengthscale of the variations of the bottom depth. When the latter are periodic, or more 
generally when they are given by a stationary ergodic process, the techniques of homogcnization theory 
are used to obtain effective long wave model equations. The two most important examples are of periodic 
bottom topography, and of topography given by a stationary random process. Recently there has been a 
renewal of interest in this problem, both from the point of view of modeling of water waves in asymptotic 
scaling regimes, and of mathematical analysis. A central question is the validity of the homogcnization 
approximation, and the character of the resulting model equations. Following [27] . the paper of Nachbin 
& S0ma [52] studies the deformation of surface waves by the effects of propagation over a rough bottom, 
taken in the shallow water scaling regime. In this work the bottom is given by a random process, and 
the authors treat both the two- and three-dimensional cases. The paper of Craig et al [10] considers large 
periodic bottom variations, again for dimensions n = 1 + d (d = 1,2), deriving model equations to quite 
high order of accuracy for the profiles which describe weak limits of surface waves in the homogcnization 
limit of the nonlinear long wave regime. Similarly, the paper of Garnier, Kraenkcl & Nachbin |12j studies 
the long wave scaling limits of water waves over a periodic bottom (for d = 1), deriving an effective KdV 
equation, for which they describe the dependence of the coefficients of nonlinearity and dispersion on the 
topography of the bottom. This study continues in Garnier, Grajales & Nachbin |13) in the case of random 
bathymetry. There are other studies of surface wave propagation over periodic bathymetry, that focus on 
regimes which are not homogenization theoretic. Namely, there is the case in which the typical wavelength 
of surface waves is comparable or smaller than the typical bottom variations. Among these, Choi & Milewski 
[6] consider periodic solutions of systems of KdV equations which are coupled through resonant interactions 
with a periodic bottom. The paper by Nahoulima et al [23] considers shallow water theory with and without 
dispersive corrections, for a periodic and piecewise constant bottom of very long wavelength. 

The paper of Grataloup & Mci |T3] considers the propagation of modulational solutions over a random 
seabed in dimension d = 1, which is extended to the case d = 2 in Pihl, Mei & Hancock [26 a . In this work, 
the typical wavelengths represented in the surface and the topography arc comparable, and the effort is to 
derive envelope equations for the free surface and to understand its statistical properties, given the ensemble 
of realizations of the random bathymetry. 

There is also a long history of study of resonant interaction between water surface waves with periodic 
bottom. The paper of Mei [5D] gives the theory of linear Bragg resonances between surface waves and 
bottom variations of the same spatial scale. This is extended to nonlinear resonances in Liu & Yue |19| . The 
difference between these references and our work is that, in the latter, short scales perturbations of the free 
surface are generated by interaction of the bottom with long waves on the free surface, a feature typical of 
homogenization theory. 

None of the references above, however, give a mathematical theorem which justifies on a rigorous basis 
the model equations that are derived. After the derivation of the shallow water model in the present paper, 
the second main point of our work is to provide a rigorous justification of this derivation. There is a history 
of results on the mathematical verification of the model equations for free surface water waves, starting in 
fact with the papers of Ovsjannikov [24] [25] and Kano & Nishida [17] which give existence theorems for the 
full water wave equations and as well a proof of convergence of solutions in the shallow water scaling limit. 
In both cases the bottom is assumed to be flat, and the authors work with initial data given in spaces of 
analytic functions. Results on long wave scaling limits of the water waves problem in dispersive regimes 
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include Craig [8] and Schneider & Wayne |28j and their treatment of the two-dimensional problem, a long- 
time existence theory, and the Boussinesq and KdV limits in Sobolcv spaces. More recently, the paper of 
Lannes |18j gives an existence theory for solutions of the water wave problem for fluid domains with smooth 
variable bathymetry, and the further paper of Alvarez-Samaniego & Lannes [2] gives rigorous results on a 
number of long wave scaling limits of the same problem (see also Iguchi [15] and earlier papers of Bona et al [4] 
and Chazel[5]), all papers working with Sobolev space initial data. In the context of this body of work, what 
distinguishes the present paper is the oscillatory nature of the bottom boundary of the fluid domain, which 
has the implications that the solutions themselves are oscillatory, and principally, that the homogenization 
theory Ansatz giving the form of solutions must be justified. Our analysis has several features in common 
with the results of |llj on the justification of the nonlinear Schrodinger equation and the Davey - Stcwartson 
system as envelope equations for modulation theory, the most important of which being that the principal 
theorem is a consistency result rather than a full fledged limit theorem for solutions. Nonetheless, as far as 
we know this is the first rigorous result which justifies with a rigorous analytic argument the application of 
homogenization theory to the water wave problem with rapidly varying periodic bathymetry. In the present 
framework, precise error estimates are needed because the shallow water limit and the homogenization limit 
do not commute. More precisely, shallow water expansions are derived for slowly varying bottoms, neglecting 
some terms that are relevant for rough bottoms. Conversely, homogenization limits are usually performed 
with low regularity estimates on solutions, that place them outside of the regime of high order shallow 
water asymptotics (see for instance [7] for a recent homogenization result at leading order for the Dirichlct- 
Ncumann operator). The point of our work and the source of many of its technical difficulties is that we 
perform the homogenization and shallow water limit simultaneously, thereby retaining the full complement of 
relevant terms from the original water waves equations. The (local) effects of this infinity of terms neglected 
in previous studies add up to create the nonlocal effects present in our approximation. 

1.1. General setting. The time-dependent fluid domain consists of the fluid domain J7(6, £) = {(x, z) <E 
R d+1 , —H + b(x) < z < ((x,t)} in which the fluid velocity is represented by the gradient of a velocity 
potential <!>. The dependent variable £(x,t) denotes the surface elevation and b(x) denotes the variation of 
the bottom of the fluid domain from its mean value. We use the Hamiltonian formulation due to Zakharov |30j 
and Craig & Sulem [9] in the form of a coupled system for the surface elevation ( and the trace of the velocity 
potential at the surface ip = &\ z=( , namely 



Writing the equations of evolution in terms of nondimensional variables, different asymptotic regimes of 
this problem are identified by scaling regimes of the associated dimensionless parameters. Denote by A 
the typical amplitude of surface waves, with A their typical wavelength. Similarly let B denote the typical 
amplitude of the variations of the bottom from its mean value Hq, with £ their typical wavelength. From 
these quantities we define the dimensionless variables as follows: 



(1.1) 




The quantity G[£, b]- is the Dirichlct-Neumann operator, defined by 




where $ is the solution of the elliptic boundary value problem 



(1.3) 




(1.4) 



x = XX' , z = H z' , t = -Aj=t' , 



( = A£, $ = ^X^TTo^ , b = Bb'(f). 
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Stemming from this change of variables there are four dimcnsionless parameters, 

n o H ° A r B 1 

(L5) M= A^' £= i/o' P = W > 7= A' 

Our analysis is concerned with the shallow water regime (i < 1. The relative amplitude of solutions is 
governed by e. In addition to this, the relative amplitude of the bathymetry is given by (3, the parameter 
7 determines the relative length of bottom perturbations with respect to the typical wavelength of surface 
waves, and the bottom variations &'(■) arc assumed to be 27r-pcriodic in all variables. We consider relatively 
large amplitude surface waves, meaning that no smallness assumption is made on e. As usual for this regime, 
we therefore set e = 1 for the sake of simplicity. With regard to the bottom variations, we set 

(1-6) /3 = ^ = 7 «1. 

The fact that (3 = 7 corresponds to small bathymetry slope in this regime, while the roughness strength is 
p := v//I/7 = 1. For clarity of notation we drop this 'prime' notation for the remainder of the paper. 

1.2. Presentation of results. The first result of this paper is the construction of an approximate solution 
(Ca,"0a) °f the water waves problem in the form of the Ansatz 

(1-7) Ca = Co(A,t) + 7 Ci(A,A/ 7 ,t/7) 

(1.8) ^ = Vo(A,t)+ 7 2 V'i(A,A/7,t/7) . 

Remark 1.1. The factor of j 2 in front of the corrector ipiis natural; indeed, this yields a 0(7) corrector 
for the velocity, which is the physical relevant quantity. 

Setting Vq = V-0o and ho = 1 + Co, we show that (Co, Vo) satisfies the classical shallow water system with 
flat bottom, 



(1.9) 



dtCo + V • (hoV ) = , 
d t Vo + VCo + (Vo ■ V)Vb = 



while the corrector terms (Ci,V'i) satisfy a linear nonlocal coupled system of equations in the fast variables 
(T = th,Y = X/i) 



{ d T Ci + Vo ■ VyCi - \Dy\ tanh(/i \D Y \)^i 
(1.10) < = Vo ■ Vysech(/i \D Y \)b , 

[ drljjx + V a ■ Vy Vl + Cl = . 

In system (jl.lOp . the functions Ci,^i are periodic in the variables Y, while the variables (t,X) are to be 
treated as parameters. The above system represents the linearized water wave equations in a fluid region of 
depth ho, with a background flow given by the velocity field Vq- The source term of the RHS is due to the 
effect of scattering of the background flow from the variable bottom. 

The second result of this paper is a mathematical justification of the derivation of the above system of 
model equations (|1.9[) (|1.10[) . Our proof is in the form of a consistency analysis of the Euler equations of free 
surface water waves, for which we show that the functions (Ca, V'a) whose constituents satisfy (|1.9|)(|1.10[) arc 
approximate solutions of the Euler equations. They are not in general an exact solution, but they satisfy 
the equations (|1.1[) up to an error term E a , and we show that this error is small. Namely, we prove that 

\E a \ H . <C 7 3/8 , 

where the appropriate norm | • \ H * is defined as \E a \ H , = \E a i\ L 2 + ~/~ 3/8 \E a2 \ H i/2, and E a = (E al ,E a2 ). 
In particular, the error is small for the usual Hamiltonian norm of the water waves equations. The most 
striking point of our analysis is that this result is valid for the natural time scale t = 0(1) associated to (|1.9p 
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only if the free surface does not resonate with the rapidly varying bottom. Such a resonance is obtained if 
there exist (t, X) such that 

(fc • V {X,t)) 2 = \k\taah(h {X,t)\k\) 

for some k € Z corresponding to a nonzero mode of the Fourier decomposition of the bottom paramctrization 
b. This condition can be viewed as a nonlinear generalization of the classical Bragg resonance which is 
obtained when the wavelengths of the free surface and of the bottom are of the same order, while here, the 
latter is much smaller. In absence of such resonances, it is possible to find locally stationary solutions for 
the corrector terms, that is, solutions to (|1.10j) that do not depend on the fast time variable r. When such 
resonances occur, the dependence of the correctors on r cannot be removed, and this induces secular growth 
effects that destroy the accuracy of the approximation (it is only valid on a much smaller time scale, t = o(l), 
than the relevant one). It is likely that in this case, the dynamics of the leading term (£o, Vq) is affected, but 
this point is left for a future study. 

The Ansatz (|1.7|) (jl.8|) and the error estimates for the quantity E a represent a problem in homogenization 
theory. The principal terms (Co^o) are solutions of an effective equation, and the multiscale terms (£1, ipi) 
are the first corrector terms. The dynamics of the Euler equations require solving an elliptic equation at 
each instant of time, on an unknown domain f2(6 7 , £) whose boundaries are defined by oscillatory functions. 
The approach we take in this paper to the analysis of this elliptic problem and its asymptotic behavior is 
to transform this domain to a reference domain Qq, resulting in an elliptic problem with rapidly varying 
periodic coefficients. The principal (effective) term and the correctors are derived from this problem, with 
the principal term solving an effective equation, and the corrector solving an appropriate cell problem. These 
are then used to express the Dirichlet - Neumann operator on the free surface of the fluid domain, which in 
turn is used to express the evolution equations . The dynamics of the short spatial scales are separated 
from the evolution of the long scales using the concept of convergence on two scales [TJ . The principal part of 
our mathematical analysis is to control the error estimates of the homogenization approximation fj 1 . Tj) (( 1 . 8|) 
in describing the solutions of this elliptic boundary value problem and the associated expression for the 
Dirichlet - Neumann problem. 

2. Euler Equations 

Zakharov showed that the water wave problem can be written in the Hamiltonian form |30j 



(2.1) dt 



C \ / I\( 5 C H 
i> ) \ -I )\ S^H 



where the canonical variables are the surface elevation £ and the trace of the velocity potential on the free 
surface if) = ^\ z=( , and the Hamiltonian H is given by 

(2.2) H(C, 1>) = \l 4>G[C, b]i/> + g( 2 dX . 

The system for (OVO is written as which in dimensionlcss form becomes 

d t C-^G^[CPb^ = , 

(2 ' 3) ^ n /JL 1 ^/!* (^[C,/^ + vC-v^) 2 

where 6 7 (-) = b{-/^) and where G^[£, /36 7 ] is the nondimcnsionlized Dirichlet - Neumann operator defined 

by 

(2.4) G^C, (3b^ = y/1 + |VC| 2 a„$ u ^ 
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and where $ is the potential function, satisfying 

f M A$ + = in n, 

in the fluid domain f2(£> 7 , £), 

n(6 7 , C) = {(Jf, ^) G -1 + 0b(X/i) < z < C(X)} . 

The operator d n is the outwards conormal derivative associated with the operator /1 A + d: 2 . One can rewrite 
in Hamiltonian form (|2.ip . replacing the Hamiltonian H given by (|2.2p by its nondimensional form 



(2.6) H (C,i') = l[ {^-G^CPb^ + C^dX . 

2.1. Notation. We denote by d = 1 or 2 the horizontal dimension of the fluid domain, and by X £ M. d the 
horizontal variables, while z is the vertical variable. We denote by e z the unit upward vertical vector. 

The domain and the potential function will depend upon both regular and rapidly oscillating variables, 
which we denote X G M d and Y £ T d = R d /(2irZ) d , respectively. That is, we will give data for the water 
wave problem which is of a multiscale nature, with the fixed multiscale bottom variations as well, and we will 
seek solutions which have a well defined asymptotic expansion in terms of multiscale quantities. To express 
this, we use the classical notation of a multiscale function that is, a function f(X, Y) defined on M. d x T d , 
for which the realization is the trace f\ Y —x = f(X,X/'y), [3]. In the problem we consider, there are other 

variables as well, such as the vertical variable z € [—1,0], for which / = f(X,Y,z) is a multiscale function 
whose realization is f(X, X/7, z). 

The differential operators V and A act on functions of the horizontal variable X. The operator A is 
defined by A := (1 - A) 1 / 2 . We use the standard notation for Fourier multipliers, namely D = AV and 

f(D)u(k) = f(k)u{k). When applied to multiscale functions, we distinguish this fact using the notation Vy, 
Ay, Dy, when differential operators act specifically on the fast variables Y, and Vx, Ax, Dx when they 
act on the long scale X variables. Finally, the notation V M stands for V = (^/J1V T , d z ) T . 

We encounter functions defined on the fluid domain f2(& 7 ,£) or the reference domain fio = K d x (0, 1), 
as well as functions defined on the free surface, parametrized by I £ R d . The notation used for function 
space norms is that j| • j| L 2, || • ||#r is used for the classical Sobolev space norms over while for norms 
defined over the boundary X £ R d we use the notation | ■ \ L 2, | • \ jjr. Norms of multiscale functions arc given 
similarly, for example | • |l 2 (c^,)- 

For all r x ,r 2 > 0, we also define the space H ri < r2 = H ri < r2 (M. d x T d ) by 

(2.7) H ri ' r %R d x T d ) = {/ € L 2 (R d x T d ), \f\ H ^.^ < 00}, 
with \f\U, r2 = |(1 - Ax)^/ 2 (1 - A y )^/ 2 /|| 2(RdxTd) . 

2.2. Change of variables and domain. The first component of the Hamiltonian (|2.6j) corresponds to the 
nondimcnsionalized kinetic energy. It follows from the definition of G^[C,(3b-y] and Green's identity that 

(2.8) / i/t-GniCPbytydX = - [ \\7^<S>\ 2 dzdX , 



where $ is the velocity potential (|2.5[) . Since this expression depends on (3 and 7 through the domain of 
integration f2(6 7 , £), it is convenient to transform it into an integral over a fixed domain independent of the 
parameters and of the perturbations £ and b. Under the assumption that the fluid height h = 1 + £ — /?6 7 is 
always non-negative, namely 

(2.9) 3a > 0, 1 + C - /% 7 > a on R d , 
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an explicit diffeomorphism S mapping the flat strip fio onto the fluid domain fl is given by 

(2 - 10) S[ {X%) Z (X,z + a(X,z)), 

where a(X, z) = (z + l)^(X) — zj3b 7 (X). We have in particular h = 1 + d z a. 

Defining <j> on Q,q by <p = $ o S, one can check (see Prop. 2.7 of [T5] and §2.2 of [2]) that the new potential 
function 6 solves 



(2.11) 



V • P[ct]V^ = 0, 



where d n <f> \ z= _ 1 is the outward conormal derivative in the new variables 

and where the matrix P[cr] is given by 

/ hi -VMVcr \ 

(2.12) i*M= T l + M|Vcr| 2 , with h=l + C~PK. 

v — h — J 

3. Multiple scale asymptotic expansions 

3.1. Ansatz and decomposition of the solutions. This section is devoted to the study of the elliptic 
problem (|2.1ip where (£, tp) are given; the time is fixed and appears as a parameter. We pose the multiple- 
scale Ansatz on (£, ip) : 

mi ( = ( (X)+rt 1 (X,X/ 7 ) 

1 ' ^ = MX) + i 2 Mx 1 x/ 1 ). 

Recalling that j3 = 7 = y/Ji, this leads to the decomposition of the height function of the fluid domain 

h = ho + (3h\ , where ho = 1 + Co an d hi = (1 — 6 7 . 
Similarly, the new vertical deformations are posed in terms of this Ansatz 

(7 = Co + /3<7i , where <jq = (z + l)(b and 01 = (z + l)(i — zb 1 . 
The coefficients P[cr] are then written as 

P[a] = P + pP u 

with 

(3.2) P = P[a ] and pP r = P[a] - P ■ 

Explicitly 



/ (Ci - 6) J -v^v^i 

We accordingly decompose the potential function <f> as 

(3.3) 4> = MX, z) + I3 1X {X, z- 7) = MX, z) + fi X (X, z; 7) 

where all the contributions coming from the roughness are contained in x- This section is devoted to deriving 
asymptotic expansions, with accompanying error estimates on the two components <f>o and %, in the limit 
/i — > 0. In order to do so, we must augment (|2.9|) with the assumption that 

(3.4) 3a Q > 0, 1 + Co > a on M. d . 

This ensures that the water depth does not vanish for the averaged fluid domain that arises when all the 
fluctuations due to the roughness are neglected. Assumption (|3 .4[) ensures the coercivity of Pq. 
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Proposition 3.1. Let (,b £ W 1 '°°(W i ) and assume that (|2.9|) and (|3.4|) are satisfied. Then for all ip such 
that Vip G H 1 / 2 (R d ) d , there exists a unique solution cj> to (|2.11|) such that V M <G H 1 (Hlo) d+1 . Moreover, O 
and x sofoe 

' V • PoV^o = 0, 
</>o |„ =0 = "00, -e z • PoV^o i^.j = 0, 



(3.5) 
and 
(3.6) 



• P[cr]V''x = -^V • Pi Wo, 
X U =0 = V<i, -e. ■ P[a]V X u = _, = ^e z • iW0o U = ^- 

Proof. The existence of a unique solution </> such that V<^ G iJ 1 (ilo) d+1 to (|2 . 1 1 1) is a classical result, and 
we thus omit the proof. Similarly, there exists a unique solution <po such that V^o G -ff 1 (^o) d+1 to (|3.5[) 
since the boundary condition on the lower boundary is the conormal derivative associated to the elliptic 
operator V M ■ PoV M . It remains to prove that \ solves (|3.6[) . A calculation gives that 

= -/3 7 e z ■ P[a]V"x U= _, - e z • P V> U= _ x - 0e z ■ i\V"0o 1^ ■ 
Since by assumption one also has 9 n 0L = _i = and — e z • PoV M 0o L = _i = 0, one has 

-e z • P[a]V"x U= _, = ^e z • PiV M <fo U^- 
It is straightforward to check that x i = ipi and that 

7P 

= -iv'-PiV^o, 

7 

and the result follows. □ 

3.2. Asymptotic analysis with estimates of V M 0o- In this section we prove an estimate on V^o, and 
we give the first terms of its asymptotic expansion in the limit as pi — > 0. For purposes of understanding the 
JT _1 / 2 -norm of the trace of V^o on the free surface {z = 0},we use L 2 estimates on fio of both V M 0o and 
its generalized Riesz transform, given by A c^V^o- This is generalized to higher order norms. 

Proposition 3.2. Let r G N and Co € W 1+r '°° n W 2 >°°(R d ) and assume that tfQj) is satisfied for some 
ao > 0. Then 

(i.) For all /i G (0, 1) and all ipo such that Vi/'o G H r (M. d ) d , the solution O to p.5[) satisfies 

WVMlz < vW — ,|Co|wi^.~ )|Wo| H r , 
ao 

HA^V^oH^ <fiC(— .IColwn+r.- ,|Cok^)|V^o|^ ■ 
ao 

(ii.) J/V^o G H r+1 (R d ) d , one also has 

||A r (V^ - V^o)|U= < ,\Co\w^ Wo\ H r+^ ; 

ao 

(iii.) Suppose that Co G W 2 + r >°°(]R d ) and A^ G ff 2+r (K d ), and set 

2 

0o 1} =- /l o(y + ^)A^o , {ho :=1 + Co). 
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as the next term of the asymptotic expansion. Then there are estimates of the remainder, in the form 
||A r (V0o- VC^o-^))!!^ <M 2 C( — ,|Co|^+--)|A^o|^ ; 

CtQ 

HA^^V^o - V-'o - M0o^)IU 2 < V 2 C( — , \Co\w-+r,oo)\Ai> \ H 2 + r . 

ao 

Proof. The first inequality of the proposition is obtained by standard elliptic estimates for || V m </>o||l2 (see 
Corollary 2.2 of [2]). The Riesz transform of V^</>o has components A^cW^ = (y^A^VS^o, A _1 df (f> ); 
estimates of the first component come from the first inequality of (i), since A -1 V is L? bounded. Estimates 
of the second component are obtained through the equation (|3.5j) itself. Namely one has an expression for 
d^(f>o in the form 



ho 



'l+/i|V(To| 2 



<9 z |Vct 



2 



(3.7) difo = u , " l2 - ' ' a^o + ^(Vao-V^o)+V-(Va ^o)-V-(/ioV^o) 



ho 



In order to get an estimate on ||A r 1 c?^ 0o 1 1 2 , we need the following lemma. 

Lemma 3.3. Let r G N and F and G > be such that A'^F g L 2 (fl ) and G G L°°((-l, 0); W^' 1 ^ 00 ^)) . 
Then 

l|A r - 1 ( I |^)IU» < C(l|G|| ir< - 1 |,o=)||A'- 1 F|U 2 . 
Proof. Just write = F — Fj%=,. Recalling that one has l/ffl/f-i < |/|ijr-i Ifflwi 7 — 1 l>°° i we S e t 

Ia-'^JIl. < IIA-^lu.ti + ll^ll^ _) 



< C(||G||,.„, l ,-„,.)||A > '-'F|| i , 



Applying this lemma to Q3.7P with G = /i|Vero| 2 , one easily gets 

||A'- 1 9^o|U= < lCo| 1 ^ + i,~)-U|A r Wo|U=, 
a y/fJ, 



□ 



and the estimate follows from the control on || A^V^^o Hi 2 established above 
Fc 

(3.8) 



For the second point of the proposition, we write <fio = ipo + /UXo '■ The resulting system for Xo ^ s 



A calculation shows that 



--V • PqV^o = -V • (/ioV^o) + a z (Va • VVo) = -/i A^o 



and that e z • PqV^Vo L=_i = 0, we obtain 



(3.9) 



xP\ m=o =0, -e z -P V^ 1) |_ 1 =0. 



Multiplying the equation by xb^ an d integrating by parts, we get 



/ V*^ • PoVxP dzdX = f Xo'^oA^o dzdX . 
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Using the coercivity of the matrix Pq (Prop. 2.3.iii of [2]), the Cauchy-Schwarz inequality and Poincarc 
inequality (in order to control ||Xo IU 2 by ||V^Xo IU 2 )' one gets 

||V"x^lU»<C(— ,|C)ki.-)|AVto|x» 

and the result follows. Higher order estimates are handled similarly and only require the control of additional 
commutator estimates. We omit these classical details. The estimate of the generalized Riesz transform 
A c^V^o — ipo) is similar to the analog estimate in (i) of this proposition. 



For the third point of the proposition, we solve 
or equivalcntly 



at lowest order in fx. We write Xq 



(i) 



with 



0o = ipo + M0o + M 2 Xo 2) 

(2) 



/io(%" + 2 )Ai/V The correction Xq satisfies the system 

V • PoV^o 2) = -V • (hoV<f> { 1} - Wod z 4 1} ) 



(3.10) 



-e z ■ P V^x|, 2) |, = _ a = 0, 



where we used the fact that e z • PoV M </> l = _ 1 = to obtain the bottom boundary conditions. Proceeding 
as above, we get the result. □ 

3.3. Asymptotic analysis with estimates of x- To find an asymptotic expansion of x, our starting point 
is the equation (|3.6p for \. Decompose the solution as the sum of a multiscale function and a correction 
term, 

(3.11) X = 4> { i\x,Y,z)\ Y ^ + y qi X [ 1 \x,z; 1 ) . 

When acting on a multiscale function of the variables (X, X/'y), the operator V M becomes Vy.z +Vj, where 

V M = ( V^Vx' 



X 



V 



V^(f(X,Y)\ Y=2L ) = [(Vy, + V» x )f{X, Y)] \ Y= , 



We can therefore write 



V ■ P[a]Wx 



so that (|3.6j) becomes (recall that j3 = ^/Ji) 



(Vy, 2 + V£)-P (Vy, 2 +V 



:(0) 



• P[tr] Vxi 



(3.12) 



V F , 2 • PoVy,.^ 05 !^ + V 7 ^ • P[a]V^xi 



(i) 



1, 

1 



+VMXi 1} )u = o = ^i^ , 

- e z ■ P Vy. 2 ^ 0) |y=f ) u= _, - VMe z • PMV^xi 1 ' 
= -e z • PiV M o | z= _! - \/^ e z ' ^U=-i 
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(o) _ p / V*^ 0) 



A = -PiW^ - P 

\ Y = 2L 

and 



In order to make the leading order terms in ()3.12j) explicit, we further decompose P as P — Pq +y/JIPfj 

h I \ n) ( -Vcr 



and Pi = P 1 (0) + ^/JIP^\ with 



and 



p (0) _ Ao/ \ (1) _ / -Va 



where the (2, 2)-coefficient of P^ is 

„(!)_,, -1/2 /Vi ^ + MlV^I 2 6-Ci\ 



Lemma 3.4. The coefficient matrix Pi = Pi(c) /ias multiscale functions as coefficients. Considered as 
Pi = Pi(X,X/y/jI), the following estimates hold for all r 6 N: 

HP^IU-h^.- + I|3 2 Pi (0) ||l~iv^ < ^ r/2 C(-, |Co|^ +1 )(Ki|^ + \V Y (l\ C r + \b\ C r + l), 

ao 

WPPhrW." + \\d Z P^\\ L? Wr^ < ^ r/2 C(—, -, \( \ C r + l, \Cl\c-, \Vx(l\cr, \V Y (l\ C r , |6| C -+0 ■ 

ao a 

Proof. The proof follows by inspecting the elements of Pi, P} and Pj . □ 

Given the above decompositions of Pq and Pi , the first term of the LHS of (|3.12l) is 

Vy,z • P Vy 2 <^ 0) = V Y , Z ■ P (0) Vy,^J 0) + V^Vy, 2 • P (1) Vy, z <^ 0) 
and, using that 7 = ^//Z, the first term of the RHS of (|3.12j) is 

I V ^.PiV^0o = iv At -Pi (o) V' i 7/> o + iv AI -Pi (o) V^(^ o ~^o) + V AI -Pi (1) V At 0o 



7 7 7 



(o) /V^oV , yfI p (o) fVV>o 

v ji^=f + v v-n ^ o 

+ -V^ • Pi (O) V^(0 o - V^o) + V • Pi (1) V^^ 

7 

Extracting the principal term from these two expressions, we deduce that 



-^The operator V* always acts on multiscale functions on the two variables X and 2 (and not on Y). The notation V M </>j 
is therefore a shortcut for V M (</>^ \y— 
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(3.13) 



Vr, 2 • P (0) Vy,^i u V = 2i + • P[<t]W X [ 



(.Ji) 



-Vy, z ■ F 
(1) 



(o) /V?/'o 







+ v^Xi )|. =0 = V>ii = . 



-e. 



(Po (0) V 



p 



(o) /V^ \ 



v^ez ■ A}|__ = _ 1 , 



with 

and with 
(3.14) 



A = A- 



7y^ 1 



V(^o - V'o) - -^PfV'Vo + (2V(t • Vy<^ 0) - ^I^^9 2 0j o) )e, 



-|Va p 



g - Vd z a ■ Vy^' - (VxCO ' V^o - (Ci - &)A^ ■ 



In order to solve (|3.13[) . we construct (p^ as a solution of a ceZZ problem in the variables F and z (the 
variable X being considered a parameter). The resulting solution cancels the higher order terms in (|3.13l) . 
and we are left with an equation for the corrector Xi ■ 

3.3.1. The cell problem. We assume that 6, £i are periodic with respect to the variable Y and we seek a 
periodic function 4>f\-,Y, z) that solves 



Vy z • P {0) V Y ,z<Pi 



(3.15) 



(o) 



h |»=o 



-01 



"Vy, z • Pi 



(0) 



d°) 



(o) /VV'o 







This choice of (f>i cancels the highest order terms in (|3.13j) . Taking into account the definition of P (u> and 



(0) 



Pi°\ we can further simplify f|3 . 1 5[) into 



(3.16) 



(h 2 A Y + d\ 

J°) 



d a ) 



= 
1 



ho 



V Y b- V^o- 



We recall that the spaces H ri ' T2 that appear in the statement below are defined in (|2.7I) . 

Proposition 3.5. The solution (j)^ of the cell problem (|3.16[) is given in operator notation by the expression 

/ q1 « ,(0) /y v \ COsh(/t (z + l)|£>y|) sinh(fe z|£>y|) Vy 

(3.17) ^ (x,y,z)= cosh(Wy|) Mx,Y)+ cosHhQlDYl) ^^■v^o(i) ■ 

Assume that ho = ho(X) = 1 + Co satisfies the hypotheses (|3.4[) . and /ei ro > d/2 and r s N. Then, for all 

multiindex a = (a 1 , a 2 ) <E N d x N d suc/i i/iai la 1 ! + |a 2 | = r, one has 

(3.18) 

ff |a,l|,r + |c,2| + |VxV , lIif°.'*0+"- + | & | C + 2 I V </>0 I r+ 1 ) , 

Moreover, derivatives of the multiscale function (f)^ | x are controlled as follows 

(3.19) ||A r V^ 0) ||i a < /i- r / 2 C(|fto|c^) x (m^IVa-V-iIh-o + |Vy^| ff o,, 0+ . + \b\ C r + 2\Vi> \ Hr+1 ), 
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and if r > 1, the same upper bound holds for -^=||A r 1 d ' z \/ ' p '<f>i \\ l 2 ■ 

It is of note that the solution of the cell problem is a multiscale expression that can be differentiated 
arbitrarily many times with respect to the variables (X, Y, z) without developing singular behavior in the 
limit as /i-)0. 

Proof. Decomposing the function <j)^ in Fourier modes tp^ with respect to the Y variable, we find 

^ = A k e hoZ ^ +B k e- hoZ ^ 

with coefficients 

1 ,. k 

e 



After substitution, the solution (f>^ is written using operator notation as in the statement (|3.17|) of the 
proposition. 

For the proof of the estimates of the derivatives of <pf (X, Y, z), the expression (|3.17p is conveniently written 
in operator notation as 

d(ho, z, D Y )MX, Y) + C 2 (h ,z, D Y )b(Y) ■ V^o(X) , 

where the components are 

r fh n x cosh(fr (z + l)|J y |) smh(h z\D Y \) Vy 

L-i(n , z, l>y) = , , , (^2[n>o, z, Uy) = 



cosh(/i |£>y|) ' v ' cosh(/i |£»r|) |Vy| 

For the first term, we remark that the Sobolcv embedding H Y ° C Ly yields 

\\9x9Y%*Ci(ho,z,DY)Ml>L- Z I [supCK^^VfeCi^i)^,-^)!^)] 2 ^, (fc = 0,l) 

x Y '" JR d z 

where V stands for Vx and V x for Vy, 2 . 

Now, by Planchercl formula (with respect to V), one easily checks that 

supKdgdfVxC^iX,-^)]^ < C(\ho\ H r+i)^ \d^V x MX,-)\^ 0+ r-wi+\^i(X,-)\ H? ^), 

1 ^ / n Y 

Plugging these inequalities into the integral above then yields the desired result, 

||^VVy,,Ci(ft 0j Z J I> y )^ 1 || £ =, x . < C(|ft | O r + i)(|VyVllHl-M.n,+l- a l + I Vy^i | ff o,. 0+ .) . 
For the control of the derivatives of C2(/io> z, Dy)b(Y) ■ V?/>0j we easily get that 

Wdidy'Vx^Mh^z^Dy^Y)- V^oll^i- < C(|/i | c , +1 )|fo|c^|V^ |^+i, 



where we have (somewhat non-optimally) estimated the action of singular integral operators on Ly at the 

cost of one derivative. This ends the proof of (|3.18j) . 

For the proof of (|3.19|) . we remark that the l.h.s. can be controlled as 
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where the summation is over (a 1 , a 2 ) G N d x N d such that | a 1 ! + | a 2 1 < r. The result follows therefore from 
(|3. 18|) and a straightforward interpolation between the lowest and highest terms in terms of /i. In order to 
prove that the same bound holds for -^=||A r - 1 <9 2 V^ 0) || L 2 when r > 1, we just have to remark that 

\\A r - l d z ^(^ 1y= ,)\\ L 2 < ^ll^ 1 ^ a (Vy I , + V^)S^ 0) || £ a £?-J 



where the summation is over (a 1 , a 2 ) e N d x N d such that la 1 ! + |a 2 | < i — 1. From the explicit expression 
of 4>\ \ we can check it is possible to replace d z (j>^ by \Dy \4>f^ in the above summation, so that the result 
follows as for ([513)1 . □ 

3.3.2. Estimate on the corrector Xi . With tf)^ as in the previous section, the system (|3.13[) reduces to the 
following boundary value problem for Xi\ 



(3.20) 



Proposition 3.6. Let r € N and denote (r — 1)+ = max{r — 1, 0} and r = (r — 1)+ + 1. The solution x± 
of (|3.20[) satisfies the estimates 

||A r V' 1 Xi 1) |U= < fi-^M r (\V7p \ Hr+ i +M r/2 |Vv,y^i|^+i,-o + \^x,y^i\h^), 

||A r - 1 a,Vxi 1) IU' < /i" 1 ^M r (|V^o|^+ 1 +M^ 1 ^ ± |Vx > yV'i|^+^a +|Vx,r^i| ff ^o+^) , 
ura'tfi M r = C(i, ^, |Co|c+inc=, ICilc+inc^ \b\c^nc)- 

Proof. The method follows the recipe of classical energy estimates, paying attention to the rapidly oscillat- 
ing coefficients and their commutators with differential operators. Indeed multiplying (|3.20[) by Xi and 
integrating by parts yields 

[ P[o-]V^Xi } ■ V^x^ dzdX = { A- V^Xi^ dzdX - f gxi } dzdX . 

J Oq J Oq J Oq 

The matrix of coefficients P[a] is coercive under condition (|2.9[) (Proposition 2.3(iii) of [2]), and is uniformly 
so with regard to the small parameters, as the scaling regime we are studying imposes that /3 = 7. 
Using the Cauchy-Schwarz inequality and Poincare inequality, one finds, as in the proof of Proposition [ 
that 

(3.21) WV^h^MoiWAWv + Wgh*) 

with Mo as in the statement of the proposition. 

For the general case r € N, the procedure is exactly the same replacing g by A r g and A by A^ r \ with 

A (r) = A r A - [A r , P[<r]] Vxi^; 
in particular, it follows from classical commutator estimates and Lemma 13.41 that 



\\A^\\ L , < \\A r A\\ L2 + C(— ,i,|Co| C ^jCi|c^J&l^+OE/ i " fe/2 ll Ar " fev ^i 1) ll^ 
so that (|3.21j) yields in this configuration 

r 

(3.22) HA^xi^ll^ < Mo(||AM|| i2 +\\A r g\\ L *) + M^^W^X^Wl^ ■ 

k=l 

We therefore need the following lemma. 
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Lemma 3.7. The following estimate holds, 

||AM|| i2 +||A' 5 || 2 < fi- r / 2 M r (\V^ \ H r +1 +^ (|VxV'l|ff r ' r o+|Vy?/>i| H r+i,r ) + |V X'>Pl\H°' r o+ r + \^Ylpl\H 1 ' r o+'-) 
Proof. - Control of \\A r A\\ L 2. Recall that 

A = -P^^ - -Pf V(<£o - V»o) - ^PfV^o 
M v A 1 

(3.23) -P ( Vx f ) + (2Va„ ■ Vy*< 0) ^l^l d J^), e , 



, ' ^ u v ' ri v ^ ho zri '\*=< 

\y = 2L 



A direct application of the chain rule gives 

||AM|| L2 < [l|Pi|U r v^.~l|A r - fc V^ 0) |U 2 + ||P 1 (0) || i « W r*.«||-A'- fc V'*(^o-^ 



L 2 

+ ||P 1 (1) |Uoo^,oc||^A r - fc V^ ||L 2 



(3.24) +||P ||^w— ||A''Vx^ 0) IIl 2 +C(|Co|c.-+ 1 ,— )||A r Vy,^i 0) || i 

"0 



Using (|2.12[) and Lemma 13.41 to control the norms of Pq and Pi , and Proposition 13.21 to control the second 
and third terms in the above expression, we find 

r 

\\A r A\\ L2 <A/ r ( M - r / 2 |VV-o|^ + i+X;M- fc/2 ||A r - fc V' 1 ^ 0) || 2 + ||A'-(Vx^ 0) | i )l|3 + l|A r (Vy )i! ^ 0) | J|| 2 ). 

We now control || A r - fc V^ 0) || 2 through (l3~19l) . while || A r (V X (^ 0) | ^ )||a and || A r (Vy, z< ^ 0) ly=x )|| 2 can 
be controlcd using (|3. 18|) and proceeding as in the proof of f|3.19D . This yields 

\Vxlpl\H r - r a + \^ X,Ylpl\H a - r o + *-) ■ 

- Control of \\A r g\\ L 2. We hrst recall that 

g = - ( V Q X ) • P (Vy z ^ 0) ) Uf + ( - Vd z a ■ Vr^ 0) - VxCi ■ V^o - (Ci - b)A^ ) i ^ . 
We get therefore 

||A''. 9 || 2 < M - r / 2 M r (|V^o|^ +1 +M r/2 ||A r (Vy, z ^ 0) 1^^)11^) 

□ 

Using (|3.22[) and the lemma, we get directly 

r 

+M r ^/i('- fe )/ 2 ||A'- fc V^ 1) || L2 , 
fc=i 

and the estimate on || A r V M Xi 1 ' ) II follows from a straightforward induction. We now turn to prove the estimate 
on ||A r_1 9 z V' t x^||. As for the control of A r_1 9*V0o m Proposition l3.2[ it is enough to get an upper bound 
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on || A r 1 9^xi. 1 ^ll- The idea is to proceed as in Proposition 13.21 using the equation to get d^Xi hi terms of 
quantities under control, 

dlx ^ = tttW V^-^t^ 9 '^ + dz{Va ■ v ^ 1)} + v ■ <y" d *xP) v • (w x «)) + RHsa 

the presence of the fast scale X/^/JI makes things a little more complicated than in the proof of Proposition 
and we need product estimates and a refinement of Lemma 13.31 for multiscale functions. 



Lemma 3.8. Let r G N, and denote (r - 1)+ = max{r - 1,0}. Let also G = G(X,Y,z) G L?w£ Y 1)+ ' 
and F G L 2 (fi ) be such that A^'^+F G L 2 (fl Q ). Then 

(r-i)+ 

W-\G ly X F)\\ L2 < \\G\\ , oow (r-i) + ,« £ ^ fc/2 !|A (r - 1)+ - fe ^|| L2 . 

k=0 

If moreover G > 0, then 

(r-l)+ 

W^YTG IIl2 " £ M-^IIA^-^+^FH^. 



Proof. The product estimates are a straightforward consequence of the chain rule; the second pair of estimates 
is derived as in Lemma 13.31 using these product estimates. □ 

Using Lemma I3~51 and the above expression for d^Xi , we g e t> with r = (r — 1)+ + 1, 
(3.25) ||A'- 1 ^Xi 1) ||L 2 < M r £ ^(v^llA^'V^lb + IIA^- 1 ^-^^^^). 

A;=0 



We can now use Lemma 13.71 to get 

\\A( r -V+- k RHS@Mh* < VTA\^ k M^ + ||A( r - 1 )+- fe 5 || i2 + ||A( r - 1 )+- fe ^|| i2 
(3.26) < /i- £ ^ i M r (|Wo|^+i +^\V x ,Yipi\ H ^ + \Vx,Yipi\H°,r+r ) + ||A (r - 1) +- fc a,A|| L 2, 

so that the only thing we still need to prove is a control on \\A < ^ r ~ 1 ' >+ ~ k d z A\\^2. 

Lemma 3.9. For all m E N, < m < r — 1. one has 

\\K m d z A\\ L 2 < fT? M r (\V1>o\ H r+i +M m/2 |Vx,v^i| ffm +i-o + \V x> yHh°^+^). 

Proof. From the explicit expression of A provided by (|3.23|) . and proceeding as in the proof of Lemma l3.7[ 
we get 

m 

||A m 9 z A|| i2 < M r ( M -^|V^ |^ +1 + ||A' m (Vx,Y^ 1 ° ) | y ^)||^L 2 +E^ fc/2 ||A m "^ z V' i ^ 0) || ff i i2 ), 
and the result follows from Proposition 13.51 □ 
The desired control on || A r ~ 1 d^Xi \\l 2 is then a direct consequence of (|3.25|) . Q3.26P and the lemma. □ 



WATER WAVES OVER A ROUGH BOTTOM IN THE SHALLOW WATER REGIME 



IT 



3.4. Asymptotic expansion of the Dirichlet-Neumann operator with estimates. We study here 
the asymptotic behavior of the Dirichlet-Neumann operator, 

(3.27) G[CPb J }^ = e z -P[a]V^ u=0 , 

where <f> solves (|2.11[) . In the previous section, we have shown that when the surface parametrization £ and 
the trace of the potential at the surface i/j are of the form (|3.ip , one can decompose <fi into 

<t> = 0o + MX) 

where </> is deduced from (|2.11[) by neglecting all the contributions due to the roughness. We have further 
decomposed 4>q and the residual \ (which contains all the roughness effects) as 

0o = V>o + ^ + M 2 x[, 2) and X = ^°) +A /7J X « 

with controls on the residuals Xo^ and Xi given in Propositions 13.21 and 13.61 respectively. We can therefore 
rewrite </> as the sum of an effective part and a residual part, 

<P = 4>eff+ A* <?Ves 

with 

<t>eff = i'o + + ^i 0) ) and 4>res = Xi X) + V^Xo ] ■ 

Similarly, we decompose the Dirichlet-Neumann operator into an effective and residual part as follows: 

Proposition 3.10. We separate the expression for the Dirichlet-Neumann operator as 

G[C,/36 7 ]^ = (G^) eff +^ 2 (G^) res 

where 

-(&l>)eff = ~ V ' (^oW'o) ~ VyCi \ Y=2L ■ VV^o 

+|£V|tanh(fto|Dy|)Vi(A:,V)| x - V^o • V Y (sech(h \D Y )b) lY=J , . 

7 7 

The remainder (Giji) res satisfies the estimate (r integer), 

(3.28) \(Gip) res \ H r < A/ rA t- r / 2 - 1 / 8 (|V^o|^+3 +^ /2 \^x.yHh^o + \^x,yHh^o+'-) 



Remark 3.11. The nonlocal operators in the expression for (Gi/>)e// arise from the simultaneous homoge- 
nization process and shallow water limit. Homogenization analysis on a shallow water expansion would give 
a different result. The reason for this difference is that certain terms neglected in standard shallow water 
expansions are not negligible in the presence of rapidly varying bathymery; their effects are described in the 
nonlocal terms of (Gip) e f / • 

Proof. Recalling that P[a] = P (0) + ^{P^ + Pi), with P 1 = pf> + y/JlP^ , one has 

G[C,/36 7 ]V = e z • P ( °V^ e# |„ + ^e z • (P (1) + Pi)V^ effU = 
+^ 3/2 e z ■ P[o-}Wcf> res u=0 
= e z • P (0) V^ e# |„ + ^e z ■ (P (1) + P^V^o |„ + ^e z ■ V"tfo U=0 
+^ 2 e z ■ (P (1) + ftJV^ + 4 0) ) |.=o + M 3/2 e z • PHV^, U=0 . 

We now decompose 



G[C,/36 7 ]V = (GV) e# +M 3/2 (G^) 
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with 

(Gi>) eff = e z • Po (0) V> e#u=0 + tfie. ■ (F (1) + P^VVo \ a=0 - 

(3.29) (G^) res = -L ez .P 1 (1) V^o|„ n +e z -(P (1) +P 1 )V^(^ 1) +^ 0) )|_ 

+e z ■ P[cr] W res u=0 . 

The two tasks of this proposition are to give an expression for {Gip) e g and to prove an estimate for (Gijj) r . 

- Explicit computation of (Gil)) e g. From the definition of 4> e ff, we have 

T^eff = ■ (^o + v^Pi (0) )V^o U=0 + e z ■ P (0) V^W u=0 + e z • P <°>W< 0) u=0 

(3.30) = -V-(/ lo V^o)-V y Ci|y = - •VVo + T^5 2 ^ 0) | 2=0 . 
Computing the last term in the RHS with the help of Proposition ^. 5[ we get 

^d^l^o = \D Y \ t&nh(ho\D Y \)MX, Y )\ Y= x + Wo ■ Vy (seeh(ho\D Y \)b)\ , 
so that G e jj is indeed given by the expression stated in the proposition. 

- Control of the residual (G?p) res . From the explicit expression (|3.29p of (Gtp) res and Proposition 13.21 and 
13.51 and Lemma HOI we get 

\{Gi>) res \ H r < ^ r l 2 M r (\VMH^+n r,2 \Vx^i\ 

r 

+M r ^ k/2 \A r ' k V^ res u=0 \ L 2 , 

fe=0 

which motivates the following lemma. 
Lemma 3.12. 

|V^ re;sU=0 | L 2 < //-^MidWol^ +/i 1/2 |Vx,yVi|// 2 -o + |Vx,yiAi|ffi.i+-o) 
Proof. We write 

|W resU=0 | L 2 < ^ 1/8 |V' i ^re.|^ o |ffi/»+M~ 1/8 |V' , 0re.|, =o |H-i/ a 

< ^V'iAV^-IU 2 + m" 1/ *I|s.v , ^.„|| £ +M" 1/8 (llv' i ^ es |U2 + ||A- 1 a,v"^ Pe a|| ia ) 

where, for the second inequality, we have used two different version of the trace lemma, namely, |Pl_ U 2 ^ 
M 1 /4|| A i/2 jF j| i2 + ^-i/^A-i/s^i?!!^ and \F lz=0 \ L 2 < \\A l / 2 F\\ L 2 + \\A-V 2 d z F\\ L2 . The estimate' follows 
therefore from the definition of 4> res and Proposition ^. 61 □ 



Bound for |A r fe V M res i |£2 are obtained in the same manner, giving rise to a power of fi in the form 

At -l/8-(r-fc)/2_ n 

Later in the consistency analysis, we will need an estimate of | (Gtp) re s Ih 1 / 2 ■ F° r this purpose, interpolating 
between H r ~ 1 and H r , we have (r > 1): 

\{G^) res \ Hr _i < A* 1/4 |(G^)r»k'+A*" 1/4 K^)re.|flr-i 

(3.31) < M P /i- (p -* )/2 - 1/8 (|V^|^+3 +^ r/2 1 Vx.yV'i l^+i-o + |Vx,y</>iUv +^) 
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4. HOMOGENIZATION WITH ESTIMATES OF THE EQUATIONS FOR WATER WAVES 

Up to this point the canonical variables of the water wave problem (C, -0) have been treated as data for 
an elliptic partial differential equation in a fixed domain. We now return to study the dynamics of the water 
waves system (|2.3[) as a time dependent problem. We first derive the effective PDEs satisfied by (Coj0o) 
and (Ci, 0i ), after which we investigate the precision of our approximate solution (|3.1[) in satisfying the full 
water wave equation. 

4.1. Effective equations. Consider the decomposition (|1.7[ 1 (|1.8[) to be an Ansatz for the full Euler equa- 
tions (|2.3[) , for which the slow and fast scale variables are identified through Y = Xj^ffx and r = t/y/Jl. 
Substituting these expressions into the full equations and using the rigorous expansion (|3.30[) . for the first 
component of equations of ()2.3|) we obtain an equation for the quantity (Ca, ipa) = (Co + v^Cu ^° + A t V'i)> 

<9tCo + <9 T Ci + ^oAV'o + VCo • VV'o + VKi ■ V0> o 
(4.1) H^v|tanh(/io|£Y|)VM^)| v _* - V0 O • Vy (sech(/i |£>y l))&| y=i = -yWi + VaW-)™. • 

At this point the fast and slow time scales are identified, Y = X/ ' y/JL and r = t/y/jl, and we have made no 
approximation. Similarly, for the second equation of (|2.3j) we obtain 

dti>o + v^r^i + Mi + Co + v^Ci + l^i'o + ^(Vy + y/Ji Vx)4'i\ 2 = 



(4.2) fi 



2 

±((G^a) e ff+ fl 3/2 (Gi> a ) res ) + (VCo + VyCl + v^ V *Cl) ' (Wo + V^Vy01 + /UVxV'l) 



2(l + /i|VCo + (Vy + ^ Vx)Ci| 2 ) 
Isolating the error terms in (|4.2p onto the RHS, one obtains 

1, 



dtlpO + y/fldrtp! + Co + VMCl + ^ | V-00 1 2 + y/W^O ■ Vy^l 

1 



2 



(4.3) = -fxdtipi - fiVr/io • ^ x4>i - ^|VyVi + v 7 ^ V^Vi| 2 

(t(( G ^)e#+ /i 3/2 (G?Aa) res ) + (VCo + VyCl + y/WVxCl) ■ (V^O + y/T^Y^l + ^ X^l)) 

+M 2(l + Ai |VCo + (Vy + 0JV A ')Ci| 2 ) ' 

The bottom profile b is function of the fast variables Y = X/ ^JJi. 

Now adopt the point of view that we seek multi-scale approximations to the system of equations (|4.1I) (14.31) . 
To impose this scaling regime, make the assumption that the variables t,X and t,Y are independent, so 
that (Co + VMCi) V'o + A'V'i) are multiscale functions of the variables (t, X, r, F). In these equations, Co and 
0o are functions of X, t only, while Ci and ip\ are multi-scale functions of both time and space. The original 
variables will be re- imposed when we return to the identification Y — X/^JJi and r = tjJJi. In order to 
justify this otherwise formal separation of slow and fast scales, we will use results on scale separation that 
appear in [3] [10]. 

Equations (|4.1[)(|4.2j) are two equations for the four unknown quantities (Co, Cij "00, "01 )■ ln order to obtain 
well defined evolution equations for them, one must identify dynamics that take place on the slow and fast 
space and time scales. This is performed using the following scale separation lemmas, in which the underlying 
periodic nature of the cell problem plays a role. 

Proposition 4.1. Let g be a continuous function on M. d which is periodic over T d ; and denote g = 
/ 2 „\d Jjd g(Y)dY its average value on T d . For any function f(X) in the Schwarz space S(R d ), we have 

(4.4) / g{Xh)f{X)dX = g f f(X)dX + 0(-y N ) , 
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for any N . 

Proposition 4.2. Let g(X, Y) be a continuous function on R d x M. d which is periodic inYG T d , and denote 
g(X) = j^jj Jjd g(X,Y)dY its average value over T d . For any function f(X) in the Schwarz space S(SL d ), 
we have 

.,N\ 



(4.5) / g(X, Xh)f{X)dX = J g{X)f{X)dX + 0( 7 iv ) , 
for any N . 

The formal derivation of the effective equations satisfied by (CojV'o) is to write (|4.1j)(|4.2[) in the sense 
of distributions, using test functions / that depend only on the large scale variable X, averaging over the 
variables Y, and neglecting all the terms that are understood to be of lower order (a rigorous justification 
of this process will be the object of the next section). Denote the mean value over the Y variables by 
(Ca^a) = (jhj J r Ca(X,Y)dY,^ J r MX,Y)dY). Then 

(4.6) dtl a = C> t C0 + d T Zl + VT^CT , 

and similarly 

(4.7) dti> a = d t 1p0 + y/Jidrlp! + fJtdtlpl ■ 

We assume that Ci = and tpi = 0, an assumption that will be shown to be consistent with equations (|4.1ip 
and (|4.12j) derived below. We obtain therefore ( a ~ Co and tp a ~ tpo at lowest order. For equation (14.11) . 
using that the K-derivative of a periodic function of Y has mean value zero, we find that 

(4.8) d t Co = -h Ai>o - VCo ■ V^o ■ 
^From equation (|4.2|) . we find 

(4.9) ^o + (Co + v^Ci) + ^|V^o| 2 = . 

The lowest order approximation to the system (|4.8[) (|4.9p takes the form of the classical shallow water system 
(with Vq = VV'o and ho = 1 + Co)- namely 



(4.10) 



dtCo = -h V ■ V - VCo ■ V 
dtV + VCo + V ■ Wo = 



Returning to (|4. 1|) ()4.2|) and using the equations (|4.10j) satisfied by (Co^o), we obtain at the next order of 
approximation the equation 

(4.11) <9 T Ci + V ■ VyCi = \Dy\ tanh{h \D Y \)tpi + V ■ V Y sech(h Q \D Y \)b . 

If Ci = at time r = 0, it remains so for all times. For the evolution equation for ipi, we find 

(4.12) dri'x + V ■ Vy Vl + Cl = • 

Again, if rpi = at time r = 0, it remains so for all times. 

The result is that (|4.10[) is the shallow water system (|1.9[) for (Co, Vo), with Vq — Wtpo, and the dispersive 
corrections are given by (|4.11[)(|4.12[) . This derivation, and a rigorous justification of it, are the principal 
subject of this paper. 
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4.2. Regularity of the approximate solution. The approximate solution is constructed from the solution 
of a system of simpler model equations (|4.10j) and (|4.11|)(|4.12|) where the first is a version of the classical 
shallow water equations. The second is the system for linear water waves in the rapid variables (Y,t), with 
a forcing term due to the presence of bottom variations, whose coefficients depend upon (Co (A X), tpo(t, X)), 
for which the slow variables are considered as being fixed. 

Theorem 4.3. For r > d/2+1, given initial data (Co(-, 0), V (-, 0)) in H r (R d ) x H r (R d ) d such that ((31} is 
satisfied. Then there exists T > and a smooth solution (Co, V ) in C([-T, T};H r (R d ) x H r {M. d ) d ) to (|4T0]) 
with this initial data. 

Proof. The shallow water system can be written as a symmetric hyperbolic system for the vector func- 
tion (Co,Vb = Vf/'o)- For Sobolev index r > d/2 + 1, these equations are locally well posed in time for 
(Co(-,0), Vb(-,0)) e H r x H r satisfying the condition (jOJ) (see for example [16]). □ 

The components of the corrector (Ci,"0i) = (Ci( t j Y\ ^ X), tpi(r, Y; t, X)) arc multiple scale functions, 
satisfying a system of the form 

(4.13) d T ( Cl )+V {t,X)-Vy ^ >( ° H^ltanh^o^X)!^!) ( Ci \ ( f 



>; W / V v.'i j \ g 

for which the large scale variables (t,X) enter as parameters. In the case of equations (|4.11j)(|4.12[) the 
inhomogeneous forcing functions are given by 

(4.14) f(Y;t,X) = V (t,X)-V Y sech(h (t,X)\D Y \)b(Y), g = 0, 

(in particular, it is autonomous, namely independent of r, and its zero Fourier mode / vanishes). 

Initial data for this system is given in the form (Ci(0, • ; t, X), ^(0, • ; t, X)) in H^(T d ) x H Y +1/2 (T d ) with 

zero mean on T d . The dependence of these solutions on the variables (t, X) will be quantified in a paragraph 

below. 

For the sake of clarity, we omit the dependence on the variables (t, X) in the statement below, since they 
only act as parameters in (|4.13[) . It is also convenient to introduce the energy or order r (r e I) defined for 
all couple of function (u, v) on T d with zero mean by 

(4.15) ||(«,«)|||r = ]T(1 + k 2 Y(\u k \ 2 + \k\t a nh(h Q (t,X)\k\)\v k \ 2 ) , 

fc^O 

where Uk and Vk stand for the Fourier components of u and v. 

Theorem 4.4. Letr e R. Forall(f,g) G C(R : H Y xH Y +1/2 (T d )) with zero mean and all (Ci(0, -),4>i(0, •)) € 

H Y (T d ) x H Y +1 ^ 2 (T d ) with zero mean, there exists a unique solution ((i,ipi) of ()4.13|) in C(R : H Y (T d ) x 

H Y +1 ^ 2 (T d )) with initial values given by (Ci(0, •), "0i(O, •)). Moreover, this solution has zero mean and one 
has 

Vrel, ll(Ci(r),Vi(T))|||.<||(Ci(0),Vi(0))|||.+r sup ||(/(r'),ff(r')|||, 

0<r'<r 



This theorem implies that (CijV'i) i s bounded in H r x H r+1 / 2 over any time interval r € [— T\,T\\. 
However, this bound may grow as T\ — > oo due to the possible presence of secular terms. Furthermore, 
when considering the dependence of this solution on the parameters (t,X), secular growth of the quantities 
dx(Ci, V'l)) d t (Ci, 4>i) is quite possible, and would affect the validity of the solution decomposition (|1.7I)(|1.8[) 
over long time intervals r E [—T/j,T/-f]. In Theorem 14.51 we show that such effects do not occur, at least 
in the absence of Bragg resonances, and for initial data (Ci(0, ■), ipi{0, ■)) chosen to be stationary in the local 
environment defined by (CoiV'o)- 
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Proof. For the sake of simplicity, we take g = in the proof below; the adaptation to the general case is 
straightforward. Considered independently of the large scale variables (t,X), the system of equations (|4.13D 
over Y £ T d has constant coefficients, and the solution operator can be conveniently expressed in Fourier 
transform. The evolution of the individual Fourier modes is described by the system 



(4.16) d T ' 



Cik \, ik-v o \ Ci* W o -u% \f Cik \ = ( h 
W v o ik ■ Vb M J W J Ui* J V o 



1 /2 

where the local velocity is Vq(£,X), and the frequency is given by w fe = {\k\tanh(ho(t,X)\k\)) ' , both of 
which depend paramctrically upon the long scale spatial variable (t,X). 

Defining new coordinates u k := cj^ 1 Ci,k an d v k := Ld^^tpi^, the propagator for (|4. 13j) (|4. 16)) is given by 



exp 



e 



ik-Vor ( cos(o; fc T) sin(w fc r) 



— sin^^r) cos(wfc7") 



-ik-Vo \ / uo k 
-ik ■ Vb / T v -Wfc 
Using complex notation for this system, define 

Z k :=u k + iv k , W k := u k - iv k , 
with which we express the general solution to (|4.16p ; 

(4.17) Zfc(r) = e- lT[uk+k - Va] Z k {0)+ f e- l{T - s ^ k+k - v ^ul 1/2 f k (s)ds 

Jo 

W k {r) = e+^-^Wj^O) + f e +t << T - s ^ k - k - Vo ^- 1/2 f k {s)ds 

Jo 

Standard use of the Plancherel identity implies that 

\\Z(T, -)\\ 2 H r + ^ < \\Z(Q, 011^+1/4 + C |r|||/( S , Olli-d-r.r]:^) 

\\W(t,-)\\ 2 k+1/ , < 11^(0, Oll^+v* + Co|-r|||/(s, Olll-ci-^]:^) , 
where we have used that uj k ~ (A:) 1 / 2 . Recovering our original variables 

l -(Z k + W k ) = u k = ^T 1/2 Cife , Y^k- w k ) = v k = u+ 1/2 il> lk 



the result is as stated in Theorem 14.41 with in addition a quantitative estimate on the growth in the fast 
time variable r. □ 

4.3. The Bragg resonance condition. Solutions to the linear equation (|4.13[) exist for all t£K, however 
(Ci("j t), t)) may exhibit secular growth in time; more precisely, it may grow linearly with respect to 
r. This is a concern for our model system because r = t/j is the rapid timescale, therefore over physically 
relevant time intervals of 0(1) in the slow time variable t, solutions of (|4.13[) may grow from 0(1) quantities 
to 0(1/7) quantities, thus leaving the range of validity of our assumption regarding the asymptotic scaling 
regime. This secular growth for Fourier modes (Ci,fcj V'l.fcX 7 ") is due to the presence of Bragg resonances of 
the k th Fourier mode of the corrector solution with the periodic variations of the bottom topography defined 
by b(Y) for which b k ^ 0. Note that these resonances differ from the classical Bragg resonances which are 
obtained with surface waves and bottoms of comparable wavelength [19] ; to our knowledge they had not 
been exhibited before. Such a resonance occurs at time t and in X if for some k 7^ such that b k =/= 0, 

(4.18) c^(X,i) 2 = (k- VoiX.t)) 2 , 

where V = VV>o, w k = (|fc| tanh(/i |fc|)) 2 . 

In absence of such resonances, it is quite easy to check that there is no secular growth of the first corrector: 
(Ci(t, •; t, X), tpi(r, -;t, X)) remains bounded with respect to r in the energy norm (|4. 15|) . This easily follows 
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from the fact that ((i, V>i) solves (|4.13j) with a forcing function / given by (I4.14|) which is independent of r. 
The time integral in (|4.17|) . which is at the origin of the secular growth, can then be explicitly computed, 
and it is obviously bounded in absence of Bragg resonance. 

Controlling the error corresponding to our approximation requires however some bound on the parametric 
derivatives of £1 and ipi (i.e. their derivatives with respect to t and X). These parametric derivatives also 
solve a problem of the form (|4.13[) . albeit with a different forcing term (f,g). Theorem 14.41 can therefore 
be used to give some control on the energy norm of these parametric derivatives; however the forcing term 
now depends on r and the linear secular growth in r that appears in the estimate of Theorem 14.41 cannot be 
removed as above. 

As previously explained, this secular growth is destructive for our approximation. While it cannot be avoided 
in general for the parametric derivatives of we still have some freedom to eliminate it. Indeed, the 

choice for the initial condition associated to (|4.13[) is so far completely arbitrary. It turns out that, in 
absence of Bragg resonance, there is one choice of initial data that removes the secular growth. This removal 
is quite spectacular since the corresponding solutions are independent of r (and therefore bounded together 
with all their parametric derivatives). These constant solutions are found by removing the r-derivative in 
(|4.13p ~ (|4.14p which leads to solving the following problem, 

V -Vy -\Dy\ t&nh(h \D Y \) \f Ci \ = ( V ■ V Y sech(ho\D Y \)b 
I Vq-Vy J \ i>i J \ 

In absence of Bragg resonance, this yields, for all k G Z, 

> (Vp ■ k) 2 sech(h \k\) . - , (K ) -fc)sech(fe |fc|) - 

K ' ' -(Vo-k) 2 + \k\tanh(h \k\) k ' ■ -{V ■ k) 2 + \k\ tanh{h \k\) k ' 

A quantitative measure of nonresonance with respect to a sequence {0 < B k < +oo : k £ Z d } is necessary 
for the analysis that follows; the k th Fourier modes (Ci.fc, $i fc) are nonresonant at (X,i) with respect to the 
homogenized solution ((o(-,t),i/jo{-,t)) and the bottom topography b(Y) if bk ^= and one has 

(4.20) \u k (X,t) 2 - (k ■ V (X,t)) 2 \ > -L . 

The sequence {Bk} is effectively a bound on the small divisor condition governing Bragg resonances, locally 
in (X,t). Given a sequence {B k : k £ Z d } such that B k < e n ^ 2 /S, if (|4T50) holds for all k ^ 0, a local 
stationary solution exists, and furthermore the secular growth of local solutions can be controlled locally in 
(X,t). When (|4.20[) holds uniformly in (X,t) for all bk ^ 0, it is a nonresonant situation (relative to the 
small divisor conditions {Bk}) and solutions can be controlled globally. 

Theorem 4.5. Let r £ N, r' > d/2 + r + 1, T > and (Co,V ) £ C r ([-T, T}; H r ' ~ r (R d ) 1+d ) be a solution 
of the shallow water equations (|4.10j) . such that 

3a >0, V(X,t) £ R d x [-T,T], 1 + { (X, t) > a . 

Assume also that the nonresonance condition j4-20\ ) holds with Bk < e h \ k \ / 5 (for some 5 > and < h < ao). 
Then there exists a unique locally stationary solution of (|4.13[) . which is given by l{4.19\ l. In particular, one 
has, for all < s <r' — r and all s' > 0, 

(4-21) \(l\c r (,[-T,T];H s x xH^) + I^HOlo ([-T,T];H*.xfff,') 

< C rss >(\( \c([-T,T\;C%), \''I>o\c([-T,T\;CZ))(\(o\ct([-T,T\;H' c ) + I V'O \c~ ([-T,T];H S X ) ) ■ 

Proof. Let F k := M x R d ^ R denned as 

(V r -fc) 2 sech((l + C)|fc|) 



-(V r -fc) 2 + |/fc|tanh((l + C)|fc|)' 
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so that C^ a ,k = F(Vo,Co)bk- The mapping Fk is smooth, vanishes at the origin, and is rapidly decaying 
together with all its derivatives as a consequence of the nonresonance assumption (|4.20p . It follows therefore 
from Moser's estimate that 

l&| S \Cl,k\c{[-T,T\;H» x ) < C{s, S ', | Co \c([-T,T];L%) , I V \c([-T,T];L%) ) (ICo \c([-T,T]-H x ) + \Vo\c([-T,T];H x ))\h\, 

so that the bound given in the lemma stems from Plancherel's inequality in the case r = 0. Bounds for r > 
and on ipi are obtained in the same way. □ 

The question as to how often Bragg resonances occur merits a discussion. For k ^ fixed, (|4.20[) is an 
open condition on the state parameters (Co, Vq) £ C . In two dimensions (that is, when d = 1), it is related 
to the local Froude number of the flow, defined by 

FAX t) -=^(M. 

' r {A ' t] - h (X,t) ' 
indeed, the nonresonance condition (|4.20[) can be stated equivalcntly as 

2 tanh(fe (X,t)|fc|) 
Fr {X ' t)= h (X,t)\k\ 

and in particular for supercritical flows Fr 2 (X,t) > 1, Bragg resonances are absent. However, the Froude 
number is an indication of criticality which is local in X , and because Vb £ H r , solutions can be supercritical 
only on compact sets. For subcritical flows, and for {X, t) fixed, at most one Fourier mode can be in 
resonanc (0. If 6(F) is a trigonometric polynomial, then there are a finite number of resonances. Any 2- 
resonances are separated by a region of non-resonance, and for Vb — > as X — > ±oo further resonances are 
avoided. In particular, if k max denotes the highest nonzero Fourier mode of b and k m i n the lowest one, then 
Bragg resonances are possible only if 

Fr 2 < Fr 2 < Fr 2 with Fr 2 = tan Hho{X,t)\k mm \) , taah(/t (X, t) \k max |) 

rain — — max ' Willi mm u f v I 5 maz t / v x M 7 I 

h (X,t)\k m in\ n,o{X,t)\k m ax\ 

For general 6(F) with infinitely many nonzero Fourier coefficients, any zero of velocity Vo is a point 
of accumulation of resonances and in particular, small resonant patches will appear because Vq — ¥ as 
X ±00. Their asymptotic strength is related to the large k- behavior of \b\\- 

The character of resonance for d > 2 is different. For b(Y) given by a trigonometric polynomial, resonances 
are isolated for the same reason as for d = 1. In the case of a general b(Y), there is the potential for a dense 
set of resonances in the state space (£0, Vo) and not just at Vo = 0. This can be seen through the parametric 
dependence on (Co,K)) of the resonant condition (|4.I8|) in wavenumber space. Given (Co,K))i this condition 
defines a hypcrsurface Ek in k £ Z d , which, if it passes through a lattice point k £ "L d with bk ^ 0, gives 
rises to a resonance. Even if it does not intersect a lattice point, under arbitrarily small perturbations at 
(Co,Vo), it will. Hence the set of resonant states (Co,K)) is dense. 

Nonetheless, in the measure theoretic sense, Bragg resonances are relatively rare. That is, there is a set 
of states (Co, Vb) for which (|4.20j) is satisfied for all k ^ 0, such that its complement has measure less than 
C5. Indeed, fix k 7^ 0. The gradient of the resonance condition with respect to (Co, Vo) on the curve Ek is 
non-vanishing, and is of amplitude of order (9(|fc|). Thus state variables (Co,^o) of distance (i5fe|fc|) _1 from 
Ek will satisfy the nonresonance condition (|4.20|) . The union over k ^ of (i5fc|fc|) _1 -tubular neighborhoods 
of the sets Ek consists of the 'bad' states, for which there exists at least one near resonance as in (|4.20|) . 
This union has relative measure bounded above by 

^\Wk< C6 - 

fc/0 1 1 



2 This follows from the fact that tanh(a:)/a; is strictly decaying on R+. 
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Therefore the resonant set is dense, but has relatively small measure in the space of states (£o, Vq). Moreover, 
we see that the set for which (|4.20p is satisfied has the character of a Cantor set. 

4.4. Consistency analysis. The purpose of this section is to evaluate the error that is made when approx- 
imating the solution (C>V0 of the full water wave problem by the functions fl3.1[) . where the components 
(CojV'o) and (CiiV'i) satisfy the effective system of equations (|4.9j) and (|4.11j) respectively. 
Write the full water wave problem ()2.3j) as 



U99) £i(C,V) = o, 

[ ' E 2 ((^)=Q. 



where E\ and E 2 identify with the LHS of equations (|2.3[) . We denote our construction of an approximate 
solution by E a := (£7i(£ a , ip a ), £?2(Cai ^a)), where (£ a ,ip a ) is defined in (|1.7j)(|1.8|) . For this approximate 
solution, the error is given by the expression 



(4.24)£ 2 (C OJ VcO = mV^o • Vj^i + ^|VyVi + VJ 1 VvV'il 2 - M^i 



2 



(4.23)S 1 (Ca,^a) = -^Tl{G^ a ) res - y/jidtt 
a) = fMVlfjo ■ Vx^l + ^|Vv^ 

U((G^ a ) eff + M 3 / 2 (G^ a ). res ) + (VCo + VyCi + yfi&xtx) ■ (V^o + \/MVr^i + MVxVi)) 
~ M 2(l + /i|VCo + (Vy + V^VxKil 2 ) ' 

The statement that the expression (Ca^a) is a good approximation for the equations (|2.3p is that the 
error is small, in an appropriate norm, for small (i. The Theorem 14.61 is a result of this form, implying the 
consistency of the approximate solution. We recall that the leading term (£o, V J o) of the approximation solves 
the nonlinear shallow water equations ([1.9)1 while, in absence of Bragg resonances, the correctors (£i, ipi) are 
explicitly given by (|4.19j) . 



Theorem 4.6. Under the assumptions of Theorem ^. 5\ the approximate solution (£ a , tp a ) given by expression 
([1.7p (II. 8| satisfies the following consistency estimates 



\El(ta,1)a)\L><C a lPl* , 
( ] \E 2 (Ca^a)\ H ^ < CaH 3 ' 4 , 

for the error term for the water wave system (|2.3|) . The constant C a is of the form 

C a = C(— ,|Co|c* , \Vo\ H * ,|6| £ »)- 

Remark 4.7. T/ie quantities (Coi V'tb Cii V^) are solutions of the model equations and following Theorems 
\4-3[ \4-4\ and \4-5\ are bounded along with their derivatives in terms of the initial data. The norm in which the 
error is measured is relatively weak, the reason being that we are dealing with a problem with fast oscillating 
functions. It is however a natural norm for this problem since it coincides with the norm of the energy 
functional associated to the water waves equations. 

Proof. The first component E\ satisfies 

\El(C a ,tp a )\ L 2 < ^/Jl\{G^)re S \L^y/^\dtC,l\L^ 

By Proposition ^. 101 the norm \{Gip) res \i,2 is bounded by /i -1 / 8 . This estimate involves norms of d, if>\ and 
their derivatives that can becontroled using Theorem 14.51 in terms of norms of the leading term £ and ^o- 
Thus the estimate of the first component of E a is shown to be as stated in the theorem. 

The second component of the error E 2 is given in ([4.24[) . It is made up of a complicated nonlinear 
expression. This nonlinear quantity consists of several types of terms, distinguished by whether they depend 
upon surface variables alone, or whether they depend upon the Dirichlet - Neumann operator and thus on 
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(4.29) Ml^lffV= < KM" 1/4 | oU= +M V4 | oUO < M 3/ V|l 2 +M 5/4 |VA| i2 +^\AVB\ L2 . 



the solution of an elliptic boundary value problem with oscillatory coefficients. Further terms of the RHS 
depend only upon surface variables, as products of functions of X and/or multiscalc functions. We will use 
the interpolation estimates in the form (for the sake of clarity, we omit the dependence on time here) 

(4.26) \f\ m/2 <M- 1/4 |/| L2 +M 1/4 |iV 
Products such as \Vipo(X) ■ Vy^i(X, X/^/Ji)\ H i/2 are controlled by 

(4.27) n\ V^o ■ Vy^i |jfi/2 < M 3/4 |VVo|ci|VyVi|ffi.-o+i 
Products of multiscale functions are bounded by 

(4-28) M||VyVi| 2 | ffl/ 2 <M 3/4 |VyVi|c^|VyVi|//o„ 0+1 

Other terms in the first line of the RHS of (|4.24[) are bounded similarly. We now turn to the second line of 
the RHS of (|4.24[) . It has the form fj,^ that we need to bound in if 1 / 2 norm. The denominator B satisfies 
B = 2(1 + |VCo + A*(Vy + y/I+V x)(i\ 2 ) > 2 and we can therefore write 

The numerator A contains many terms. To bound its L 2 norm, we have for example terms of the form 
(4-30) |(I GV ,) e//) 2| L2 < C 

where C depends on |V(o|c 2 i I^V'olc 2 ; \^ x,yC,i \h 3 ' t q+ 1 i \^x,y i>\ \h 3 - t o+ 1 ■ Here again, Theorem 14.51 is used 
to control the last two quantities in terms of norms of £o and ipo . Estimates of terms of the numerator A 
which depend upon the quantity (Gijj) res from the Dirichlet - Neumann operator use the results of Section 3 
on the boundary value problem with periodic oscillatory coefficients. For example, 

(4.31) |(M 1/2 (GV0r es )V <Cy/ 4 

Examination of all terms leads to \A\li < C. Noting that the computation of WB gives one factor /i and 
that each derivation costs a factor /i 1 / 2 , we find 

(4.32) |VB| L 4 < <V /4 

Considering all terms similarly, we arrive to the conclusion of Theorem 14.61 □ 
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